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Abstract. We present a posteriori error estimates for a recently developed atomistic/continuum 
^— ( coupling method, the Consistent Energy-Based QC Coupling method. The error estimate of 

the deformation gradient combines a residual estimate and an a posteriori stability analysis. 
The residual is decomposed into the residual due to the approximation of the stored energy 
and that due to the approximation of the external force, and are bounded in negative Sobolev 
^ norms. In addition, the error estimate of the total energy using the error estimate of the de- 

formation gradient is also presented. Finally, numerical experiments are provided to illustrate 
our analysis. 
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1. Introduction 



Quasicontinuum (QC) methods, or in general atomistic/continuum coupling methods, are a 
class of multiscale methods for coupling an atomistic model of a solid with a continuum model. 
These methods have been widely employed in computational nano-technology, where a fully 
Ch atomistic model will result in a prohibitive computational cost but an exact configuration is 

required in a certain region of the material. In this situation, atomistic model is applied in 
^-H the region which contains the defect core to retain certain accuracy, while continuum model is 

^ applied in the far field to reduce the computational cost. 

QQ A number of QC methods have been developed in the past decades and are classified in two 

T^J- groups: energy-based coupling methods and force-based coupling methods. Despite the fact 

that the force-based methods are easy to implement and extend to higher dimensional cases, 
CN energy-based methods have certain advantages. For example, the forces derived from an energy 

potential are conservative which could leads to a faster convergence rate in computation, and 
the energy of an atomistic system can also be a quantity of interest in real application. However, 
^ consistent energy-based coupling methods can be tedious and restrictive on the shape of the 

• ^ coupling interface in more than one dimension (see [11^ |3] ) and it was not until recent that a 

^ practical consistent energy-based coupling method was created by Shapeev |10|, which is the 

^ Consistent Energy-Based QC Coupling method that we analyze in the present paper. 

A number of literature on the rigorous analysis of different QC methods have been proposed 
since the first one by Lin f5]. However, most of the analysis are on the a priori error analysis, 
and only a few are on the a posteriori error analysis. Arndt and Luskin give a posteriori 
error estimates for the QC approximation of a Frenkel-Kontorova model [H O I]- A goal- 
oriented approach is used and error estimates on different quantity of interests, each of which 
is essentially the difference between the values of a linear functional at the atomistic solution 
and the QC solution, are proposed. The estimates are decomposed into two parts, one is used 
to correctly chose the atomistic region and another is used to optimally choose the mesh in the 
continuum region. Serge et al. [9] give error estimates, also through a goal-oriented approach, 
of the original energy-based QC approximation, whose consistency is not guaranteed. Both 
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of the above works employ the technique of deriving and solving dual problems as a result 
of the goal-oriented approach. Ortner and Siili [7] derive an a posteriori error indicator for a 
global norm through a similar approach as ours. However, the QC method analyzed there does 
not contain an approximation of the stored energy which is essentially different from the QC 
method we are interested. 

The present paper provides the a posteriori error analysis for the Consistent Energy-Based 
QC Coupling method [ lOj for a one dimensional periodic chain with nearest and next nearest 
neighbour interactions. The formulation of the QC approximation has the feature that the finite 
element nodes in the continuum region are not restricted to reside at the atomistic positions, 
which creates the situation that interaction bonds often cross the element boundaries, which 
is common in two dimensional formulation. We then derive the residual in negative Sobolev 
norms and then the a posteriori stability constant as a function of the QC solution. The error 
estimator of the deformation gradient in L2-norm is then obtained by combining these two 
analysis. In addition, we derive an error estimator for the total energy difference by using that 
of the deformation gradient. It should be remarked that though both of the error estimators 
are global quantities, they consist of contributions from element. As a result, an adaptive mesh 
refinement algorithm is developed and applied to a problem that mimics the vacancy in the 
two dimensional case, and the numerical results are presented. 

1.1. Outline. In Section we first formulate the atomistic model through both a continuous 
approach, i.e., the deformation and the displacement are considered as continuous functions on 
the reference lattice, and a discrete approach, which is always taken in previous literature. We 
then formulate the Consistent Energy-Based QC Coupling method in one dimensional setting. 

In Section [Sj we derive the residual estimates for the Consistent Energy-Based QC Coupling 
method in a negative Sobolev norm. The residual is split into two part, one is due to the 
approximation of the stored energy and the other is due to the approximation of the external 
force. 

In Section |4j we give the a posteriori stability analysis. 

In Section [5| we combine the residual estimate and the stability analysis to give the a pos- 
teriori error estimate of the deformation gradient in L2-norm and that of the total energy. 
In Section |6j we present a numerical example to complement our analysis. 

2. Model Problem and QC Approximation 

2.1. Atomistic Model. As opposed to taking only a discrete point of view in many QC 
researches, we use both continuous functions and discretized vectors to denote the displacement 
and the deformation. The reason for doing this is that the Consistent Energy-Based QC 
coupling method, which we analyze in this paper, is easily formulated through the continuous 
approach, while discrete formulations could make the residual analysis of the external forces 
much easier. 

For an infinite reference lattice with atomistic spacing e, we make the partition = 
{^/}£-oo of the domain M such that M = U^_.^r/ and T/ = [{£ - l)e,fe]. We then de- 
fine the displacement and deformation of this infinite lattice to be continuous piecewise linear 
functions u,y£Vi{T^)n C°(M). We use u and y to denote the vectorizations of u and y such 
that ui = u{ie) and yi = y(fe). We know that ui and yi are the physical displacement and 
deformation of atom £ respectively. 

To avoid technical difficulties with boundaries, we apply periodic boundary conditions. We 
rescale the problem so that there are N gN atoms in each period and e = 1 /N, which implies 
that u and y are 1-periodic functions and u and y are A^-periodic vectors. We also impose a 
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zero-mean condition to the admissible space of displacements, which is defined to be 

U = {uePi{T')nC°{R):u{x + l) = uix)sind [ u{x)dx = 0}. (2.1) 

Jo 

The set of admisible deformations is given by 

y = {yeViir')nC\R):y{x) = Fx + u{x),ueU}, (2.2) 

where F > is a given macroscopic deformation gradient. 

As we mentioned above, it is necessary in the analysis of the external forces to employ the 
discretization of the displacement and the deformation. Therefore, by the relationship between 
u, y and their vectorizations u, y, the discrete space of displacement and the admissible set of 
deformation are defined by 

N 

W = {ueR^ :ue+N = ue,e^ue = 0}, (2.3) 

i=i 

and 

ye = {y^M.^: y^^^ = Fie + u^, w G W}, (2.4) 

where the zero-mean condition on the displacements, i.e., e X^^i = is obatined by applying 
the trapezoidal rule to evaluate the integration u{x) dx with respect to the partition T'^ and 
using the periodicity of u. 

For simplicity of analysis, we adopt a pair interaction model and assume that only nearest 
neighbours and the next-nearest neighbours interact. With a slight abuse of notation, the stored 
atomistic energy (per period) of an admissible deformation is then given by 

uv) ■= e f: ^(^'^^' - - + .f: ^c'^^' - - ^'^>) 

i=i e=i 

N N 

e=i e=i 

where G C^((0, -|-oo)) is a Lennard-Jones type interaction potential. We assume that there 

exists > such that (j) is convex in (0, r*) and concave in (r*, -|-oo). 

For the formulation of the external energy, we first define the linear nodal interpolation 
operator : C°(M) ^ ViiT') n C^{R) such that 

legiis) = g{ie) V(? G C^M). (2.6) 

Then given a dead load f gU, we define the external energy (per period) caused by / to be 

„i N 

if, u), := / Ie{fu)dx = V efiue =: (/, (2.7) 
Jo 

where / and u are the vectorizations of the external force / and the displacement u according 

to 

Thus, the total energy (per period) under a deformation y G ^ is given by 

E,{y;F)=Uy)-{f,u)e, 

as u is determined by y and F. However, in our analysis, we always assume that F is given 
and as a result, we simply write Ea,{y; F) as Ea,{y). 
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The problem we wish to solve is to find 

ya G argminE^iy), (2.8) 
where argmin denotes the set of local minimizers. 

2.2. Notation of Partitions, Norms and Discrete Derivatives. Though it is natural 
to introduce the QC approximation after the atomistic model, we decide to pause here and 
introduce some important notation that are used throughout the paper in order to make the 
flow of the paper more smooth and save some space. 



In Section 2.1 , we have introduced the partition of the domain M. We now fix the notation 
for a generalized partition. 

Let T"" = {r^}^_oo be a given partition such that TJP = where > xf_^ are 

the nodes of the partition. We denote the size (or the length) of the fc'th element by := 
= ^T~^k'-i- We also define the mesh size vector e'^ such that := (£^)^_oo S (M+)^. 

Given a partition and a function g G C^{R), we define the Vi direct interpolation 
Im : C°(M) ^ Pi(r™) n C°(M) by 

iI^g)ixT) = gixT) V5GC0(M), (2.9) 

and Img is often denoted by gm- We also denote the vectorization of g £ C^(M) with respect 
to T™ by g"^ such that 

9T = 9ixT). (2.10) 
Let P be a subset of Z. For a vector v € and a partition 7""*, we define the (semi-)norms 

= ( (E,.per|H')'", l<P<oo. 
[ maxgfzv \vi\, p = oo. 

In particular, if rim is the number of the nodes of 7"™ that are in [0, 1], we simply define 



\V\\(^p 



max£=i ...„„Jr;£|, p = og. 



We now define discrete derivatives. Suppose v G C'^(M) and v"^ is its vectorization according 
to 7"™". We define the first and second order discrete derivative v"^' by 

= -^-^^ and ,v-l = (2.11) 

wheie£f:= ^{e'j + e'}_^_^). 

It can be proved that for v"^ G C''(M) n'Pi(7~'") and being its vectorization, we have the 
identity 

\\v""\\Lr'm = h""\k^- (2.12) 

Since is special and uniform, we simply use e and e to denote its mesh size vector and 
mesh size without superscripts and subscripts. 

In addition to these, we denote the left and the right limit of an open interval u by L^^ and 
Ruj, which are also used later in our analysis. 
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2.3. QC Approximation. The QC approximation we analyze in this paper is essentially the 
Consistent Atomistic/ Continuum Coupling method developed in flOj. We briefly redevelop this 
method in ID so that it is easily understood and enough for us to carry out the analysis. 

We first decompose the reference lattice, which occupies M, into an atomistic region Q^, which 
should contain any 'defects', and a continuum region f^c, where the solution is expected to be 
smooth. Moreover, we assume ^Ig, to be a union of open intervals and Jlc to be a union of closed 
intervals, and ^l^D^lc = M. Since we impose periodic boundary conditions on the displacement 
and notationally it is easier to assume the atomistic region is away from the boundary of the 
period we analyze, we make the following assumptions on J7a and 0,c: 

• and appear periodically with exactly period of 1, i.e, if x G ila then x + 1 G f^a 
and the same for il.c- 

• 36 > 2e such that ((0, 6)U {1 — 6, 1)) C Qc, i-e., the atomistic region is contained in the 
'middle' of the chain. 

Note that, there is no such a restriction, that the interfaces where different regions meet should 
lie on the positions of the atoms, i.e., it is not necessary that ilc H r^a G eZ, which was always 
assumed in previous modeling and analysis of ID QC method. 

Then, in order to reduce the number of degrees of freedom, we make the partition 7''* = 
{r^}^_oo of the domain M according to the above region decomposition of M as follows: 

• = and Tk+K = [x'^ + l,x^_-^ + 1] = [x^_^^, which implies that 
the partition is i^-periodic with | Tf^\ = 1, and there are K elements in each [0, 1]. 
We also assume that is the left most node and x'^ is the right most node in [0, 1]. 

• If £e G Oaj then 3i € Z such that Xj = ie, i.e., every position of an atom in the atomistic 
region is a node of this partition. 

• dQc is a node in this partition which means that each element is contained in only one 
of the two regions. 

• l-^fc I ~ ^fc — '^k i-e., the size of each element in the continuum region is 
larger than or equal to 2e. 

We emphasize two definitions 

4 := max{^ : £e < x^} and Ok := (2.13) 

which are extensively used in the analysis and significantly simply the notation. Note that 
0<9k<l. 

Based on this partition of the domain, the QC space of displacement and the QC set of 
admissible deformation are defined by 

Uqc = {ueVi{T'')nC°{R) :u{x + l) = u{x) and [ u{x)dx = 0}, (2.14) 

Jo 

and 

3^qc = {y e ViiT'') n CO(M) : y(x) = Fx + u{x),u£ U^,}. (2.15) 
The discrete QC space of displacement and the QC set of admissible deformation are defined 

by 

< = {u^ G : 4 = 4+K,yk G Z, and -(x^i - xti)4 = O}, (2.16) 

k=l 

and 

y^, = {y'' G : = Fxk + ulu€ U^,}. (2.17) 
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Note that unhke and , in which every vector has the physical displacements and defor- 
mations of the atoms as its components, lA^ and only contain vectors whose components 
are the values of displacements and deformations at the nodes of T^. 

The approach to couple the atomistic and continuum energy is to associate the energy with 
interaction bonds. The term bond between atoms i G Z and z + r G Z refer to the open interval 
h = {ie, {i + In our case, since only nearest neighbour and next nearest neighbour bonds 
are taken into account, r = 1, 2 only. 

To develop the coupling method, we define the operator D^y for an open interval uj = 
{Luj,Ruj) CR and y G C^{R) such that 

D^y:=-^{y{R^)-y{L^)). (2.18) 

\UJ\ 

If we take any y G C°(M) n T^i'°^(M) deformation (note that this 'deformation' might 
be non-physical) and a bond b = {ie, {i + rb)e), we can define the atomistic energy contribution 
of bond b to the stored energy to be 

a,(y) = l!lM^{rtDi,nn,y), (2.19) 

and its continuum energy contribution to the stored energy to be 

Cb{y) = - [ HVrMx))dx, (2.20) 
'"b Jbnnc 

where Vr^y = rfey'(x). 

Since we are only interested in the situation in [0, 1], which is extended periodically to the 
whole domain, the set of bonds that we will consider is 

B = {{ie, {i + r)e) : r = 1, 2, i = 0, 1, . . . , iV - l}. (2.21) 

Therefore, coupling the two energy contributions together, the stored QC energy (per period) 
of a deformation y G C''(M) H M^^'°°(M) is then given by 

£cic{y) = Y,[ab{y) + Ci,{y)\, (2.22) 

b&B 

which was shown in jtlOj to be a consistent coupling method, where the definition of consistency 
is as follows: 

£'^{Fx)\v] = 8'^^{Fx)[v] = G C°(M) n VF^'~(R). (2.23) 

Given a dead load f G U the QC approximation of the external energy (per period) caused 
by / is given by 

{f,ur,)h := / Ih{fuH)dx = .(4+1 - 4-i)ft4 =■■ {f\u'^)h, (2.24) 
•^0 k=i ^ 

where Ih is the linear nodal interpolation with respect to , and and are the vector- 
izations of / and u/j. 

Thus, the total energy (per period) of a deformation y^ G 3^qc is given by 

Eqc{yh;F) = £qc{yh) - {f,uh)h- 

For the same reason that F is given, we write Eqc{yh'i F) as Eqc{yh)- The problem we wish to 
solve is to find 

yqc G argminEqc(3^qc). (2.25) 
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3. Residual Analysis 

In this section, we bound the residual in a negative Sobolev norms. We equip the space U 
with the Sobolev norm 

ll'^llwi.2 = II^^'IIl2[o,i]' foruGZY, 
and denote it by lA^''^. The norm on the dual := (U^''^)* is defined by 

||r||^-i,2 := sup T[v], iorTGU-^'^. 

1|i'1Imi.2=i 

In the following sections, we formulate the problems in variational forms and then analyze 
the residual. 

3.1. Variational Formulation and Residual. Let ya be a solution of the atomistic problem 
(2.8). If ya!{x) > on [0, 1], £a.{y) has the variational derivative at ya and therefore, the first 



order optimality condition for ( |2.8| ) in variational form is 

£L{y.)[v] = {f,v)e V^GZY, (3.1) 

where 

'?a(ya)M = eY,(t>'{nDby^)nDt,v. (3.2) 



Let yqc be a solution of the QC problem (2.25). If yqc'(x) > on [0, 1], then <fqc(y) has the 
variational derivative at yqc and the first order optimality condition for (2.25) in variational 
form is 

£'(ic{y(ic)[vh] = {f,vh)h yvheUqc, (3.3) 

where 

£'qciyh)[vh] = [a'b{yh)[vh] + cb{yh)[vh]] 

b&B 

|6 n ^a\(t>' {nDbnn^yh)Dbnu^Vh + — / (1^' nyh)^ n^h dx Mvh S Z^qc. 

(3.4) 



In conforming finite element analysis, where the finite element solution space is a subspace 
of the original solution space, the residual is defined as the quantity we obtain by inserting the 
computed solution to the equation which the real solution satisfies. However, in our case, 3^qc is 
in general not a subspace of 3^a) and hence the functional £a,{-) is not defined on 3^qc in general 
and £qc{ ) is not defined on 3^a either. The way through which we circumvent this difficulty is 
to define mappings between the solution spaces so that the residual could be well defined. In 
concrete, we define Ju '-U ^ l^qc and Ju^^ : l^qc U such that 



-Iqc 

K 



Jwn = 4n - ^ ^(x^+1 - x^_i)u(x^) Vu G ZY, (3.5) 

1=1 

and 

TV 

^Wqc^^fe = huh - e X] ^^/^ ^ ^qc- (3.6) 



1=1 
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It is easy to check that Jyu and Ju^c'^h satisfy the corresponding mean zero condition of U and 
Uqc, which imphes that Juu £ Uqc and Ju^^Uh G ^- With a slight abuse of notation, we define 



Juy = Fx+Juu = Fx+hu-^ Y.^4+i-4~i)u{xk) = Ihy-\ ^[xl^i-xl_^)u{xl) \fy G U, 

(3.7) 



i=i e=i 
and 



1 ^ 1 ^ 

Ju^^Vh = Fx + Ju^^Uh = Fx + leUh - 2^ X] '^hi^-e) = hy - 2*^ ^ ^Vh G U^c- (3.8) 

£=1 £=1 

We then define the residual (at the solution yqc) to be 

= E'^{Ju^,yc^c)[v] - £;qc(yqc)H 

= [^a(<4/qcyqc)H - {f,v)e] - [£'q^{yc^c)[Juv] - {f,Juv)h] 

= KiJu,.ycic)[v]-£Uy^c)[Juv]] + [{f,Juv)h - {f,v)e]. (3.9) 

and understand i? as a functional in ly(~^'^. By this formulation, we essentially split the residual 
into two parts: the first part is the residual of the stored energy and the second part is the 
residual of the external force. We will bound these two parts in the following sections. 

3.2. Estimate of the Residual of the Stored Energy. In this section, we analyze the first 
part of (3.9), which is the residual of the stored energy. 

Before we give the theorem, we make several definitions that simplify our notation. 

First, we define the set fCc to be 

)Cc:= {k:ke{l,...,K} such that rfcn[0,l] / but rfcn(l,+oo) = and Tk C n^}, (3.10) 

which is essentailly the set of indices of the elements in the continuum region in [0, 1]. 

Second, suppose the atomistic region consists of M disjoint subregions in [0,1], i.e., f^a H 
[0, 1] = ufi^Ql^ among which 0!^ fl f^i = if i 7^ j, we define the nodes lie on the atomistic- 
continuum interface of the atomistic regions be x^^_, i = 1, . . . ,M and those lie on the right 
interface be ic^^., i = 1, . . . , M. 

Third, we define /C^ C /Cc to be the set of indices of the elements in the continuum region but 
not adjacent to an atomistic region, i.e., Vfe € /C^, k / La^ and /c — 1 / Ra^, Vi G {1, 2, . . . , M}. 

Using these definitions, we have the following theorem. 

Theorem 1. For yh & y with y'f^{x) > 0, we have 
where 



>2 Z^^L-r . 2 

j=0 j=0 

2 , 2 



= (^E[['^']]L,+. + rT.ii^']]l-.^^)K (3.13) 



2 ^ 

j=0 j=0 
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if k = La^ for some i £ {1, 2, . . . , M}, i.e., is adjacent to and to the left of an atomistic 
region, and 

2 2 

= (^E[['^']]L.+. + ^^E[['^']]4+.)^ (3-14) 

j=0 j=0 

if k — 1 = for some i G {1,2,..., M}, i.e., Tk is adjacent to and to the right of an atomistic 
region. [[(I)']]e 's will be defined in the proof. 



Proof. By (3.2) and (3.4), we have 

KiJu^cVhJiv] - £'qciyh)[Juv] =e (l)'{rbDbJu^^yh)rbDbV - ^ 1^ n {rbDbnn^yh)Dbnn^Ju^ 



beB ■JbnUc 



=e X] ^'i'^bDbyh)rbDbV - ^ |& n irbDbnn^yh)DbnnJhV 



beB 



beB 



^ — I (t>'{Vr^yh)VrJhvdx, 



(3.15) 



since DbJu^^Vh = Dbhyh = Dbyh, D^Jyv = D^hv for any w being an open interval, and 
{Juv)' = {Ihv)\ which can be easily verified by noting that JuqcVh and Juv are leVh and IhV 
shifted by some constants. 

To make further analysis of (3.15), we subtract and add the same terms 

|& n ^i,\(t>' {i'bDbnn^yh)Dbnn^v and / (I)' {V ri,yh)^ nV 



beB 



beB 



to get 



beB 



Ki'^Uqcyh)[v] - £qciyh)[Juv] = e(t>' {nDbyh)rbDbV - \b n n^lcp' {rbDbnn^yh)Dbnn^v 

i>'{Vrtyh)VrtVdx 



fb JbnQc 

E |6 n (rbDbnn^yh) [DbnnJhV dx - Dbnn^ 



beB 



I E - / ^'i^nyh)['^nlhv - Vr,v] dx]. (3.16) 



We first analyze the second and third groups, which turn out to be as we will see immedi- 
ately. 

For the second group, we have, 



'^\bn (rbDbnn^yh) [DbnnJhV dx - Dbnn, 
beB 

2^ |6n i2a|(?> [rbDbnn^yh) 



beB 



RbnUe, — Lbnfie, 



Rbnn^ — Lbnn^ 



(3.17) 
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We define the above to be if 6n l^a = O- If 6n Jla 7^ O, since both -Rbn^a -^^fenric ^-re either 
at atomistic postions in 0,^ or on dQc, they must be nodes in T^. Therefore, by the definition 
of IhV, the following holds 

Ihv{LbnnJ = v{Li,nnJ and Ihv{RbnnJ = v{RbnQj, 



which implies that (3.17) is 0. 

For the third group, upon defining xs to be the characteristic function of a set 5, we can 
rewrite it as 

= X] 7 X] / Xb(t>'{^r,yh) i^nhv - Vr^v] dx 
beB keJCc 

2 If, 
= XI X] X] 7 / ^b(p'C^ri,yh) ["^rJhV - Vr^v] dx 
r=l b&B,n=r k&Kc '^^ 

= Y.Y1 Yl - XbC^'iVrVh) [Vrhv - Vrv] dx 
r=l k£Kc b&Bn=r '^^'^ 

2 r 1 



(3.18) 



r=l fce/C, 



bGB,ri,=r 



since 'VrVhlTk is a constant on each element. By the ID bond density lemma flOl Lemma 3.4], 



~a.e. 1) 



bGB,ri,=r 



we have 



2 r I 

"/"'(^ryhln) [ Yl -Xb][VrIhV-Vrv] 

2 r 

= E XI ^'^^rt,yh\n) r{lhv{xk) - hvixk-i)) -r{v{xk) - v{xk-i)) 

r=l fcG/Cc 

Again by the definition of IhV, 

hv{xl) = v{xl) and Ihv{xl_^) = v{xl_^), 



(3.19) 



and thus (3.19) is 0. 

Now we turn to the analysis of the first group and analyze 

e(t)' {nDbyh)rbDbV -\hr\VLi,\ct)' {rbDbnQ,^yh)Dbr]Q,^v - — [ (j)' {V r^yh)V r^v dx (3.20) 

^b Jbnnc 

for each interaction bond b. 

If 6 C fia, we have IftnOd = rb£, l&nfial =0 and the equivalence of the operators Db = -DfenfJa- 
We know that (3.20) is by substituting these equivalences. 
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If 6 C n Tk for some k G /Cc, then |6 D 0,c\ = and \b D f^al = 0. We also note that 
^nUhix) = nDbUh, as yn is affine on Tk, and ^ /^^j^^ Vr^v = enD},v. Using these equivalences, 



we know that (3.20) is again 



Therefore, we only need to analyze the bonds crossing the atomistic-continuum interface 
or the boundaries of two adjacent elements in Og. Because of its tediousness, we leave the 
detailed analysis to the Appendix but just present the result here. Employing the notation 
often adopted by a posteriori error analysis for elliptic equations, we have the following result 

£ (j)' {rbDbIeyh)rbDbV - ^ \b D ns,\(f)' {rbDbnQ^yh)Dbnn^v + ~" / (f>' {^nyh)^nV dx 
i=l ^ 

M . 

(3.21) 



i=l 



where for k = Lai 



(3.22) 



c/>'((i - oMn^, + euy'h\n) - (i - h)<t^' {y'h\n^.) - Ok<t^' {y'h\n) 



+ 



(1- W( 



-yh\n+2 + 



2(1 -e, 



-yh\Tk+i 



) 



+ 



+ ek<l^'{2y',,\T,) - ct^'{y'h\n+, + (1 - Gk)y'h\n+, + Gky'h\n) 
(1 - 0fc)0'(2y;,|T,+J + ekck'{2y'^\T,) - <t>' {{I - ek)y'H\n^, + (1 + Ok)y'h\n) 



(3.23) 



b']]i+2 = '/•'(^Y^y^y/ilTfe+i + Y^y'h\n+2) - <l^'{y'h\n+2 + (i - ^k)y'h\n+^ + ^ky'hln), 

(3.24) 



for k = Raj 



[[<P% = 'A'((i - &k)y'h\n+^ + Gky'h\n + y'h\n.^) - (t> 
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29h 



-yhln + 



yh\n^^) + a-W{'^y'h\n^.) 



and for k £ IC' 



+ 



1 + 0^— ' 1+C7fe 

(/-'((i - Ok)yh\n+^ + Gky'h\n + y'h\n.^) 

m]i = ^'{'^yh\n) - </''((! - Ok)y'h\n^, + (i + ^fc)2/;iTj, 
(1 - 9k)4>'{y'H\n^,) + ^fc</''(y;.lTj - </''((! - efc)y;,|T,+i + e^y'^W,) 
2(1 - 0fc)0'(2y;,|T,+O + 2M'(2y;.|Tj 

<A'((2 - ^fc)y;.|T,+, + euy'h\n) - 'A'((i - Ok)y'h\n^, + (i + ^fc)yhlT,) 



(3.26) 
(3.27) 

(3.28) 



(3.29) 

= <A'(22/;,|t,+J - <A'((2 - ^fe)y;.|r,+, + Oky'^W,). (3.30) 

Distributing the contribution of ( |3.21[ ) to each element and applying Cauchy-Schwarz inequahty, 
we obtain the estimate stated in the theorem. □ 

3.3. Estimate of the Residual of the External Force. We now turn to the estimate of 
the residual of the external energy. Upon defining Jyv := Vh^ the residual of the external force 
is given by 



{f,Vh)h - {f,v) 



(3.31) 



where f,v gU. 

To further analyze (3.31), we introduce a new partition = {TJ}^^^ of the domain M, 
such that all the nodes in partition and partition are included in this partition. The 
indexing of the nodes in follow the rule that the node in is labeled as Xj^ in T^. We 
also assume there are n nodes in in [0, 1], i.e., 

n = I {e, 2e, . . . , Ne} U {xi , X2, . . . , xr} \ , 

where |^| denote the cardinality of a finite set A. 

The inner product associated with partition is then defined by 



{f,9)r 



„1 n 

/ /.(/5)dx = j;-(x,Vi 
^0 2 



(3.32) 



where Ir is the linear nodal interpolation operator with respect to T*", and and are the 
vectorizations of / and g with respect to . 

Now we decompose the residual of the external force into three parts by adding and sub- 
tracting the same terms, 

(/, Vh)h - if, V)e = [if, V)r - if, V),] + [(/, Vh)r " (/, v)r] + [(/, Vh)h " (/, Vh)r] • (3.33) 
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The following three lemma are derived to give the estimates of the three parts. 

according to and . 



Lemma 2. Let f,v,f^,v'^ be the vectorizations of f,v G 
Then the following inequality holds 



la 



1 



\{r,vnr-{f,v)e\<^e'\\f'hl^lC^)\\v'\k., 



(3.34) 



where KLu = {A; G {1, . . . , K} : Xk 7^ ^k^}, in other words, ICu is the set of indices of the nodes 
in such that x^ does not coincide with any of the nodes in T'^ . 

Proof. We first write out the two inner products and eliminate the terms that are the same 

"1 ^1 
{f,v'')r - {f,v)e = ^el-{flv1 + -^e-ifivi + fl+lVl+i) 



= e and Z^+iW^+i = flj^+iV^^+i, i = 1,2 



, ^, . . . , 



(3.35) 
4, if 



as - £jk+3 

ikE ^ Xk and 4+ie / Xk+i- 

For k such that i^E 7^ x^, by the definition of /, v, and , we have /g^ = fj^-i, 
= v^j^_^, = and = uj^+i. We als o hav e fj^ = (1 - 6'fc)/4 + 6'fc/4+i and 



Vj^ = (1 — 6k)vef. + 9kVi^,+i. Inserting these equalities, (3.35) can be estimated as 



(r,^''>r--(/,^)e| 



+ ^(1 - ek)efe,+ivi^+i + ^(1 - ^fc)4(l - ^fc)/4 + [(1 - Gk)vi, + ^fc^^^+i] 

I E ^4 [^'^(^'^ - - /4)^4+i] - M0k - l)(/4+i - /4)^'4 

^ Je3|0fc(l-^,)/4+l<+l| 



kelCu 
< 



^^(E-i/4.ii^)YE^K.ii^)'^r'ii^' 



8 , 

■ k&Ku 

which concludes the proof 



\q{Ku)\\'" Up 



(3.36) 
□ 



Remark 1. If /C = 0, i.e., every node in is also in T"^, then this part of the residual is 0. □ 

Lemma 3. Let f,v£ C°(M) n'Pi(7'^) and Vh = IhV G C° (R) OViiT^) be the Vi interpolation 
of V according to partition. Let f^,v^ and be the vectorizations of f,v,Vh respectively 
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according to , and fCc is defined in (3.10). Then we have the following estimate 

X] ^k\\f\\'ei^(vl) 



{f,VH)r-{f,v)r = {f,vl)r-{ry)r< 



k=Kc 



(3.37) 



e'i + e'i+i), hk = Ujk+i - jk)£ and Vl = {jk + 1, . . .,jk+i - !}• 



Proof. Using the fact that {vDj^ = v'j^ and by Cauchy-Schwarz inequahty, we have 

n 



ifc-1 

k&ICcj=jk-i 

^ E f E 



E 

i=ift-i+i 



(3.38) 



where = 2 (ej + ^j+i)- Upon defining gf such that gj = {v 



h)3 



(note gj^ = gj^^^ = 0) and 



by Lemma 15 in Appendix [C] (Discrete Friedrich's Inequality) and Rieze-Thorin Theorem, 



ife-i 
i=jfe-i+i 



j=jk-i 



i=ifc-i+i 



E <oafc-J^-iH E ^^^f ' (3-39) 



where g'j = ^ = {f^^Yj — i'^hyj- ^ appears in the last inequality since maxj < e. Since 

and are both piecewise linear on T*", we have (fj^)j — {v^Yj = {v' — v'f^){x) Vx G (3;J_]^,xp, 
and as a result, 



Ik 



i=ife_i+i 



|(i;' -U;,)(x)Pdj; = - t'^||L2[^r ] 

■Ik — l -Ik 



lk-1 



By Lemma 12 in Appendix [ 



X ,X ■ 



Put all the results above together and apply Cauchy-Schwarz inequality, we obtain 



< 



r ife-l 

E E ^K/J)^ 

ifc-i 1 



i=ifc-i+i 



E E ^1/; 



kGlCc 



i=ifc-i+i 



< 



jjj— 1 

E^K E ^/I'^ 

fc='Cc j=jk-l+l 



The eatimate in the theorem holds as lit;'! 



L2[o,l] 



\v'y for G C70(M) n Pi(r^). 



1^^ IIl2[o,i]- 

(3.40) 
□ 



Lemma 4. Let f,v£ C°(M) nPi(7'^) and Vh = hv G C°(M) nPi(T'') be the Vi interpolation 
of V according to the T^. Let f^,v^ and v"^ be the vectorizations f,v and Vh according to . 
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If Vh = 0, then we have the following estimate 



{f,Vh)h - {f,Vh)r < 



k=Kc 



+ 



r/||2 



k=Kc 



£2 



(3.41) 



where ICc is defined in (3.10), Vj. = {j^ + 1, . . . ,jk+i} and will he defined in the proof. 



Proof. Since Ihifvh) is also piecewise linear with respect to the partition, we apply the 
trapezoidal rule here to evaluate {f,Vh)h = Jq Ihif^h) dx to obtain 



{f,Vh)h - if, Vh)r 

= E A^.+ih{fvh){xl) + E + e]+i)h{fvH){x^j) + AWvhX^^ 



i=ifc-i+i 



1 1 

j=jfc-i+i 



•. (3.42) 



We define g and G such that = {fvh){x'j) and = {Ih{fvh)){x'j). It is easy to check 
that gj^. = Gj^ and 



G 



idjk - 9jk-i) G ICc and i = 1, . . . , jfc - jfe-i, (3.43) 



where = Yl]^=j^. 1+1 = Xk — Xk-i- Therefore, by Theorem 16, we obatin the following 
estimate 



\{f.Vh)h-{f\Irv'')r 



keKc 



1 ^ ^ 



j=ifc_i+i 



^ E 4 I 



(3.44) 



where is the second finite difference derivative with respect to the T^. 

By the definition of and vl, gj = ifvh){x-) = f^{vl)j. Using {vD'^ = Vj E Vl g] can 
be written as 



5" = (K);' = (r);'K). + |(r);-(^D;- + %i(r);-+i(^D;-+i- 



(3.45) 
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Noting that ^ < 2 and < 2 and defining hk : 
the following estimate 

(/, Vh)h 



( (jk-jk - 1 )e) (ifc - ifc - 1 +!)£) 



we have 



-1 J/C ^ 



< 



fce/c 
1 



i=jfc_i+i i=ife-i+i 

1 1 XV ff II 1 1 r 1 1 I /1 1 1 xTf 

11/ lle(o^)ll^/^lle(c?,) 



" 1 I 



T \/ 1 1 



r//||2 



+ 



E ^tiiri 



r\/| 



(3.46) 



For further estimate, we first bound ||i'^||^2^ by ||(i'/i)'||£2^. Since Vh{x) is piecewise linear 
with respect to 7^ partition, we can apply the trapezoidal rule to the integration on each 
element to get 



T.'M)j = T.2^''i+'U)(<h = T.'^J2^(''h)j + ivlh^^^ / Vh{x)dx = 0. (3.47) 
j=l j=l j=l -^0 



The last equality holds by the periodic condition on v^- Thus, we can apply Lemma 14 
Appendix [C] and Riez-Thorin Theorem to obtain 



m 



^ .3 = 1 ^ 



(3.48) 



Since v'f^{x) = {vlYj on (x^_i,xp. 



Whfdx = ||w^||l2[o,i]. 



By Lemma 11 in Appendix [B| 

ll'^hllL2[0,l] ^ lb'llL2[o,l] = \\v'\\i2. 

Combine these results, the estimate stated in the theorem is easy to establish. 



(3.49) 



(3.50) 
□ 



Having the three lemma and distribute the contribution to each element, we now give the 
theorem which essentially gives the estimate of the residual due to the external force. 



Theorem 5. For f,v(^U and Ju defined in ( |3.5[ ), we have 

f2, 1 



\\U.Ju-)h - <{Y.4Y=-- ^ext(/), 



(3.51) 
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where 

+ ^(-)^/^tiiir'ii,V(^.)+/.tiiirii,\,^^ 1', (3.52) 



and JCc is defined in (3.10), P| is defined in Lemma^ hk, is defined in Lemma^ and hk+i 
is defined in Lemma^ 

Proof. We can not directly apply the three lemma to estimate the three parts in (3.33). The 
reason is that Jkv 7^ Vh, which is the direct interpolation of v according to T^. The way to 
circumvent this difficulty is by defining 

w:=v-'^ -{xk+i - Xk-i)v{xk) and Wh := hw = Jyv, (3.53) 
k=l 

and noting that 

(/, Juv)h - if, V)e = if, Wh)h - if, W)e " (/, C), = (/, Wh)h " (/, w) „ (3.54) 

and w'{x) = v'{x) Vx G M. Then by the three lemma, we have 

\{f,Juv)h - {f,v)e\ = \{f,Wh)h - {f,w)e\ < { Yl '^fc ^ } ' H^' lU^ [0,1] = { ^fc ^ } ' II ^' lU^ [0,1] , 

kaK-c k&fCc 

(3.55) 

which establishes the estimate in the theorem. □ 

4. Stability 

Stability of the QC approximation is the second key ingredient for deriving an a posteriori 
error bounds. Since we would like to bound the error of the deformation gradient in L^-norm, 
we derive the stability estimate in this section. The procedure of deriving the a posteriori 
stability condition largely follows that of the a priori stability condition in [8]. 

For an a posteriori error analysis, the natural notion of stability for energy minimization 
problem is the coercivity(or, positivity) of the atomistic Hessian at the projected QC solution 
Ju^cVqc- 

KiJu^cyqc)[v,v] > Ca(2/qc)|b||i2[o,l] Vi; G U, (4.1) 

for some constant Ca(yqc) > 0. To avoid notational difficulty, we vectorize the above inequality 
and work on W instead. Let Ju^cV'^^ be the vectorization of JuqcVqc, then (4.1) is equivalent to 

E':iJu,^y^'')[v,v] > Ca(Jw,.2/^^)||^||| G U. (4.2) 



In the remainder of this section, we derive the explicit condition on y'^'^ such that (4.2) holds. 
The Hessian operator of the atomistic model is given by 

N N 



=1 i=\ 

|2 



We note that the 'non-local' Hessian terms + t^^^-^l can be rewritten in terms of the 'local' 
terms \v'^ and and a strain-gradient correction, 

/, / |2 r)l'|2|0| / |2 2|//|2 
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Using this formula, we can rewrite the Hessian in the form 

N N 

E':{y)[v,v] = eY,Mve\' + eY,Bi\v';\', 

e=i 1=1 

where 

A,{y) = 4J\y[) + 20"(y;_, + y[) + 24J\y\ + y^+J (4.3) 

sKy) = -</'"(y; + ym). 

Recah our assumption that is convex in (0, r*) and concave in (r=|,,+oo). For typical pair 
interaction potentials, y'^ < r*/2 can only be achieved under extreme compressive forces. Since, 
under such extreme conditions a pair potential may be an inappropriate model to employ 
anyhow, it is not too restrictive to assume that the the projected QC solution Ju^dJ^^ satisfies 

{Ju,.y''% > n/2 e z. 

As a result of this assumption, and the properties of (p, we have —(p"{y'^ + > G Z and 
thus Bi>0\/i£Z. 

As an immediate consequence we obtain the following lemma, which gives sufficient conditions 
under which the a posteriori stability of QC approximation can be guaranteed. 

Lemma 6. Let Ju^cV^^ ^ satisfies mmi{Ji(^J'f^ > r*/2; then, 

E':{Ju,^y'^'')[v,v]>A,{Ju,^y^'')\\v'\\%^ ^veU, where A^iJu.^y^'') = ^min^Ae{Ju,^y^'). 



The coefficients Ag^Jn^^y'^^) are defined in (4.3). 
Proof. If min^(JiYq<,y^^)^ > r*/2, then 



N N 



e=i 1=1 

Before we present the main theorem and its proof in the next section, we state a useful auxil- 
iary result: a local Lipschitz bound on E'l. The proof of this Lipschitz bound is straightforward 
and is therefore omitted. 

Lemma 7. Let y, z ^ such that min^ y'l^fJ- and min^ z'^ > /j, for some constant /z > 0, 
then 

where Cup = M^{[n, +oo)) + 8M3([2/i, +oo)) and Mi{S) = max^e^ W{S)\. 

5. A Posteriori Error Estimates 

5.1. The a posterior error estimates for the deformation gradient. The error we esti- 
mate is e := ya — JuqcVic the W^'^-norm, as for yi,y2 G y, yi — y2 G ^- To avoid technicalities 
associated with the nonlinearity of our models, we make an a priori assumption: we assume the 
existence of the atomistic and QC solutions and make a mild requirement on their smoothness 
and closeness (cf. (|5.1[)). 



Theorem 8. Let y^c be a solution of the QC problem (2.25) whose gradients are such that 



min^ [JUc^cy'^'^)'e — G Z and A^{Ji(^^y'^'^) > 0, where A^ is defined in the statement of 
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or 



Lemma^ Suppose, further, that is a solution of the atomistic model (2.8) such that, fi 
some T > Q, 

\\{y^ - Ju,.ync)'\\L^[o,i] = 11(2/" - Ju^.y^^yhr < ^- (S-i) 

Then, if t is sufficiently small, we have the error estimate 

\\y^ - Ju,.y^c\\mOA] = 11(2/' - Ju^.y^'l'lk =^ A^){^-torc{y^c) + ^cxt(/)), (5.2) 
where the functional of the residual of the stored energy (^store(') is defined in ( |3.11[) and the 



functional of the approximation error for the external forces ^ext(') is defined in (3.51) 



Proof. From the mean value theorem we deduce that there exists 6 £ coiw{y^, Ju^cV^'^} such 
that 

= {£Uy^)[v]-£',{Ju^^y^^)[Juv]) 
-{{f,v)e-{f,Juv)h). 



The first group was analyzed in section 3.2 Theorem [T] and the second group was analyzed in 



section 3.3 Theorem [Sj Inserting these estimates we arrive at 

EM[e,e] < (<rstorc(y^^) + <^cxt(/))||e'||,2. (5.3) 
It remains to prove a lower bound on £i[{6)[e, e]. From our assumption that mm{Ji( y'^'^)'^ > 



r=K/2, and from (5.1) it follows that 

min0« > r*/2 — r. 

e 

Assuming that r is sufficiently small, e.g., r < ri := ^ m.mf{Ji{^^y'^'^y^, we can apply Lemma [t] 
to deduce that 

£':{9)[e,e] >£:{Ju^^y^^)[e,e] - Cu4iO - yn\\i^\\e'\\% 

>S':iJu,.y''l[e,e]-Cupr\\e'\\l, (5.4) 

where CLip may depend on ri. 

We can now apply our stability analysis in Section |4| Since {JuqcV'^^ye ^ for all £, 
Lemma [6] implies that 

£:{Ju,.y'^'')[e,e]>MJu^^y^^)\\e'\\%, 



which, combined with (5.3) and (5.4) , yields 



{A.{Ju,.y''l - CLipT)||e'||| < £':{e)[e,e] < (^storc(yqc) + 4xt(/)) ||e'||,|. 

Dividing through by ||e'||^2 , and assuming that r < min(Ti, r2) where T2 = A^:{Ju^^y'^'^) / {2Cup) , 
we deduce that 

^^^%^||e'||,2 < (<fstore(yqc) + <^ext(/)), 

which concludes the proof of the a posteriori error estimate for the deformation gradient. □ 
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5.2. The a posterior error estimate for the energy. Besides the deformation gradient, the 
energy of the system is another quantity of interest. In this section, we derive an a posteriori er- 
ror estimator for the energy difference between the atomistic model and the QC approximation, 
namely, 

^a(y'')-^qc(yqc). (5.5) 

To analyze this difference, we decompose ( |5.5[ ) as 

\E^iy'') - -Eqc(yqc)| = |Sa(y^) " ^a(Jw,,yqc)| + |^a(Jw,,yqc) " ^qc(yqc)|. (5.6) 

We then analyze the two groups separately. 

To analyze the first group, we have the following Lemma. 

Lemma 9. Let y,z €z y and y,z G be their vectorizations, such that min^?/^ > fi and 
mini -2^ ^ A* for some constant /U > 0, and y G argmin£'a(3^)- Let e = y — z, then 

|i?a(Z/)-i^a(^)| <Cfip||e'|||, (5.7) 
where C^p = ^Miiii^, +oo)) + 2M3([2/i, +oo)), where M,{S) = max^e^ W{S)\- 

Proof. We first rewrite the difference of the total energy as the summation of the differences of 
the stored energy and that of the external energy: 

E,{y)-E,{z) = {8,{y) - 8,{z)) - {{f,z), - {f,y),) 

For the difference of the stored energy, we have 

N N 

^ 1 

N N 

e=i e=i 

where S,} £ convjy^, z'^} and ^| e convjy^ + y^+i, 4 + 

For the difference of the energy caused by the external forces, we have 

{f,z),-{f,y), = -{f,e), = -£'M[e], 

by the first optimality condition of y € argmin£'a(3^)- It is then easy to obtain the estimate 
stated in the Lemma by using Cauchy-Schwaz inequality to the non-local term. □ 



Lemma 10. For yh G 3^qc and y'}^{x) > 0, we have 



where 

1 



Vk 2 / J LLT jjt-fc-i-r_/ ' 2 
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if k £ K,'^, 

1 1 ^ 

^l=Y. + 2 E (5.10) 

i=-i i=-i 
if k = Lai f^''" some i £ {1, 2, . . . , M}, and 

1 - 1 



E + 1 E (5.11) 



if k = Ra^ for some i G {1, 2, . . . , M}. [[</']]£^ 's and t]^^, 's will be defined in the proof. 
Proof. We first decompose the energy difference to two parts: 

Es.{Ju^,yh) - Eqcivh) = [£a.{Juc,^yh) - £qc{yh)) - ((/, Ju^,yh)e - {fh,yh)h)- 

We first analyze the energy difference of the stored energy. Since ri^DhJu^^yh = ^bDbyhi we 
have 

£s.{Ju^cVh) = ^ab{Ju^,yh) = ^e(l3{ri,Di,Ju^^yh) = ^e(l){rbDbyh), (5.12) 
fees beB 6eB 

and 



^qc(y/i) = E ['^b{yh) + Cb{yh)] = E ^'/'('^fe^fenoay) + ^ / (t>{'^ rt,y{x)) dx 

We analyze the energy difference bond by bond, 
If 6 C r^a, then 

ab{yh) + Cb{yh) = Obiyh) = e4'{rbDbyh) = ab{Ju^,yh), 

and the energy difference in this bond is thus 0. 

If 6 C n Tfc for some k S JCc, then |6 n Oa| = and 



(5.13) 



abiJuqcVh) - [abiyh) + Cb{yh)] =ab{Ju^,yh) - Cbiyh 

Tb JbnQ 



e<p{rbDbyh) - ^ I (?!>(V,.^y/i(j;)) dx 



[(l){rbDbyh) - H^nyhix))] dx. (5.14) 



_ 1 

n Jb 

Since yh is affine on T^, Vr^y{x)) = VbDbyh and subsequently, ( |5.14 ) is 0. 

We are left with the interaction bonds crossing the atomistic-continuum interface and the 
boundaries of the elements in the continuum region. Again because of its tediousness, we leave 
the detail of this analysis to the Appendix and only give the results here: 

M 1 Ml 1 

uju^^vh) - £M = E E + E E MW.,+j + E E MW+r (5-15) 

1=1 j=^i i=i j=-i keK.'^j=-i 

For k = Lai where i S {1, 2, . . . , M}, we have 

MV, =e{<P{{l - ek)y'h\T,^, + 9ky'h\n) " (1 " OMy'^b,^,) - eM^\T,)], (5.16) 

M]t,_, = 2e{m - ek)y'h\n+, + (1 + Ok)eky'h\n) - (i - 0M2y'h\n+.) - (i + eM^y'^^W,)] ■ 

(5.17) 
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and 

- (2 - ^.)<^(^y;.iT,+. + - Okm'H\n)}. i^-m 

For A; = where z G {1, 2, . . . , M}, we have 

= e{<t>{{i - ek)y'h\n^, + hy'hW,) - ek<P{yh\n) - (i - Ok)4>{y'h\n^,)] , (5-i9) 

16 2 

- (1 - ^,)</)(2y;,|T,+j - (1 + ek)4>{^^y'^\T, + Y:p^y;.k_0 }, (5.20) 

and 

= 24 '^((2 - Ok)y'h\n+, + ^ifcy;.lTj - (2 - 0fc)</.(2y;,|T,+ J - e^4>{^y'^\T,)}. (5.21) 

For k e K'^, we have 

[['^]]4 = - Wh\n+, + eky'^ln) - (1 - ^fc)<^(y;.lT,+J - eMhln)}, (5.22) 
[[<^]]4-i = l^im^ - ek)y'h\n+, + (1 + ^fc)yklTj - (1 - ^fe)</'(22/klT,+J - (1 + ^fc)0(2y;,|Tj}, 

(5.23) 

and 

[M14+1 = ^420((2 - Ok)y'h\n+, + ^fey^lTj - 0fc^(2y;,|T,+J - (2 - 9k)cl>{2y'^\n)} ■ (5.24) 

We then analyze the energy difference caused by the external forces. The energy difference 
is given by 

if, Ju^c'^hle - if, Uh)h = if, Uh)e - if, Uh)h, 

since Ju^c'^h = Uh + C for some contant C and {f,C)e = VC. We decompose this energy 
difference to each element and write it as 

K 

\{f,Uh)e-{f,Uh)h\<^vi:k^ (5.25) 
fe=i 

where 



4fe = 



(l-^fc)^e(/4-i'"4-i + /4-i+i^4-i+i) + ^ Yl £{fm + fe+iue+i) 

^=4-1+2 

1 ^ 1 

+ 6'^£(/4^4 + /4+i^4+i)} - 9 (^fc - a^fc-i) [/(xfc_i)y(xfc_i) + /(xfc)y(xfc)] 



2' 

fe=i 

(5.26) 



where = f{is) and = n(fe). □ 
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6. Numerical Experiments 

In this section, we present numerical experiments to illustrate our analysis. Throughout this 
section we fix = 1, = 8193, and let (f) be the Morse potential 

4){r) = exp(— 2a(r — 1)) — 2exp(— a(r — 1)), 

with the parameter a = 5. 

For our benchmark problem, we defined the external force / to be 

We briefly explain the meaning of the external force. On each atom, the external force is a 
product of three components. The third component, namely — ^ is essentially ^ where 

ri is the distance between an atom and the center of this atomistic chain located at +0.5. 
This non-linear force will create a defect in the middle of the chain but affect little in the far 
field. The second component, namely 1 ]v-t — > adds a decay of the first component and in 

particular, it is when £ = N, which prevents the 'kink' of the force on the boundary due to a 
rapid change of the sign of the force that will leads to non-smooth deformation gradient that 
should be contained in the atomistic region. The first component, which is the constant 0.1, is 
to rescale the force so that the solution of this problem is stable. 

We solve for the atomistic problem and consider the solution to be the accurate solution. We 
then solve for the QC problem on different meshes generated by the mesh refinement schemes. 

We show two relative errors against the number of degrees of freedom. The first one is the 
error of the deformation gradient in L2-norm over the L2-norm of the difference between the 
deformation gradient of the atomistic solution and the homogeneous state, which is defined by 

'^deformation • n , 7-1 ii ' v^"-*^/ 

The second relative error is the absolute value of the energy difference of the atomistic solution 
and the QC solution over the absolute value of the energy change of the atomistic solution from 
the homogeneous state, which is defined by 

._ l^^aCya) -^qc(2/qc)| 2) 



\EM-E^{Fx)\- 

Before we present the plots of the errors, we first introduce the mesh generating schemes. 



6.1. Mesh Construction. To avoid unnecessary technical difficulty in the mesh refinement 
algorithm, we assume that the defect core is already captured in the middle of the chain. There 
are three mesh generating schemes we use. 

The first mesh generating scheme is derived in Section 7.1 of [6j using calculus of variations. 
From this analysis, we get that the (quasi-)optimal mesh size in the continuum region, with the 
restriction that the atomistic region is symmetric and has K atoms on each side, is given by 



(6.3) 
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Since the mesh size can not change continuously and we restrict the smallest mesh size in the 
continuum region to be 2e, we use the following algorithm to generate this mesh (we only list 
the case on the right hand side of the atomistic region) : 



Algorithm 1. 

(1) Set atom + 1 to be the middle of the atomistic region. 

(2) Choose K so that there are K atoms on each side of the atomistic region. 

(3) Choose h to be 2e for every element on the right hand side of the atomistic region until 
h{r) > 2e, where r is the distance between the right boundary of the previous element 
and the middle of the atomistic region. 



(4) Choose h according to (6.3) until the right boundary of the newly created element is 
out of the right limit of the chain. 

□ 

The second mesh generating scheme is essentially a mesh refinement process according to 
the error estimator with respect to the deformation gradient according to Lemma [2| Lemma [3] 
and Lemma [4j The mesh refinement algorithm is stated as follows: 



Algorithm 2. 

(1) Set atom + 1 to be the middle of the atomistic region. 

(2) Choose 5 atoms on each side of the atomistic region. 

(3) Divide the left and the right part of the continuum region into two equally large element 

(4) Compute the QC solution on this mesh and then compute the squared error indicator 
of each element rji and sort these indicators according to its value. 

(5) Bisect the first M sorted elements such that 

A/-1 M 

rjf < 0.5??2 and J^Vi > 0•5^?^ (6.4) 

i=l i=l 

where r]i is the error estimator of each element defined by 

^deformaUon ^ ^(^e)2 ^ (^/)2] ) i / A4^^|cj^ _ (g g) 

If the element is near the atomistic region, merge the element into the atomistic region. 

(6) If the resulting mesh reaches the maximal number of degrees of freedom, stop the 
process, else, go to Step 4. 

□ 

The third mesh generating scheme is the mesh refinement process according to the error 
estimator with respect to the energy which is defined by 

^energy ^ ^E^^^deformaUony ^ ^ ^/^^ ^g g^ 

for each element and the refinement algorithm is exactly the same. 

In short, the first and second mesh generating schemes tend to minimize the error in the 
deformation gradient and the third one tends to minimize the error in the total energy. 
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Figure 2. Efficiency Factor of the Gradient 



6.2. Numerical Results. We compare the relative errors defined in 6.1 and 6.2 We plot the 



relative errors against the number of degrees of freedom with respect to the meshes generated. 

Figure [T] shows that the pre-defined optimal mesh performs better than the two mesh refine- 
ment strategies for a fix number of degrees of freedom. The possible reason for this is that, due 
to some technical difficulty in coding, both of the mesh refinement algorithms tend to produce 
larger atomistic region by merging the elements in the continuum region to the atomistic region 
and create some unnecessary degrees of freedom. For the two mesh refinement strategies, the 
one according to the gradient error indicator perform better asymptotically. 
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Figure [2] shows the efficiency factor of the error estimator of the deformation gradient. It 
shows that the efficiency factor is comparatively large but decreases as the number of degrees 
of freedom increases and finally become stable. The reason for this phenomenon lies in the 
form of the external force. One can show that if the external force takes the form of /(r) = ^, 
where r is the distance to the centre of the defect, then the residual due to the external force 
is of order /i^ as opposed to order h in general which is achieved by our analysis. As a result, 
our estimate exaggerate the real error by ^ for this particular external force. This phenomenon 
gradually disappear as the continuum region moves apart from the centre of the defect since the 
influence of this exaggeration is eliminated as the external force tends to when it is away from 
the centre of the defect, which makes the residual of the sotred energy become the leading error 
term. It can also well explain the fact that the efficiency of the estimate is better for the mesh 
refinement strategies than the pre-defined mesh for a certain number of degrees of freedom, as 
the two mesh refinement algorithms tend to put more atoms in the atomistic region, i.e., the 
continuum region is further away from the centre of defect than that of the pre-defined mesh. 



^ 10"^ 
a 

s 
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Figure 3. Relative Error of the Total Energy 



Figure [3] shows that the refinement based on the energy error performs the best among all 
the three mesh generating schemes. 

Figure [4] shows the efficiency factor of the error estimator of the energy. For the same reason, 
this factor decreases as the number of degrees of freedom increases and finally becomes stable. 



7. Conclusion 

We have presented the a posteriori error estimates for the Consistent Energy-Based QC 
method in one dimension. The procedure of the estimate is the same as that in [8]. However, 
since the formulation of the QC problem is newly developed and is totally different from previous 
ones, new techniques have been developed and applied to deal with the difficulty in the analysis. 
Several results derived may be of independent interest and usefulness. In addition, the error 
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estimate of tlie total energy is also derived. Numerical experiments are also implemented to 
illustrate our analysis. 

Particular interesting future work are the extension and the implementation of the a posteriori 
error estimate in higher dimensional problems. The difficulty lies in the complication of the 
formulation and the varied location of the interaction bonds. However, since a priori analysis 
for the two dimensional problem has been proposed |j6j, ways of circumventing these difficulties 
could be a source of reference. 

Appendix A. Detailed Analysis for the Residuals of the Stored Energy 

In this section, we provide the omitted detailed analysis for the residuals of the stored energy, 
namely 

£L{Ju^cyh)[v] - E'^^{yh)[Ju[v\ and S^{Ju^^yh) - SqciVh), 

where yn G 3^qc, y'^ix) > Vx G M and v eU. 
The idea is to find the differences defined by 



e(l)'{rhDhyh)rbDhV - \bnn^\<l)'{rhDbr,n,yh)Dhnn,v / (/)'(Vr,y/i)Vr,u dx, (A.l) 



for each interaction bond b. 

We have analyzed the cases that 6 G and b £ Tk Ci Clc and are left with the analysis for 
the cases that b is across the atomistic-continuum interface and the boundaries of the elements 
in the continuum region. There are three cases and in each case there are three subcases to be 
considered. 

Case 1: b is across two adjacent elements Tk,Tk+i G ^c- In this case |6n J7a| = and the 
atomistic contribution of the interaction bond in the QC energy is 0. 




and 
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Subcase 1: If 6 = (4e, (4 + l)^) , then = 1, 6 n Tfe = [4£, x%h^ Tk+i = [xj^, (4 + l)^], 
nDbV = v'^^_^_^ and 



nDbVh = = (1 - OkjVhWk+i + QkyhWk- 



We have 



{nDbyh)rbDbV - — [ cj)' {V r^yh)^ r^v dx 

=e(f)'{{i - 0k)yh\Tk+^ + dkyh\Tk)ve^+i - ^(P'{nyh\n+i) [ nv'dx - ^4>'{nyh\n) I nv'dx 
=eU'{rbDbyh) - (1 - ek)4>' {y'hW,^,) - 0k4>' {0ky'h\n)\v[,+^, 



and 

£(t>{nDbyh) - — I 0(Vr(,y(a;)) dx 

=£0((i - 9k)y'h\n+, + dky'hln) - I 4>{y'h\n) dx - [ Hy'^k,^,) dx 
=e(^ci>{{i - 9k)y'h\n+, + Gky'hh) - dMln) - (1 - 0kmy'h\n+^)) ■ (A.3) 

Subcase 2: If 6 = {{i - l)kS, (4 + , then = 2, 6 n = [(4 - l)e, x'^],bn Tk+i = 
[xt {h + , nDbV = v'^^^i + v'^^ and 



nDbyh = = (1 - dk)yh\n+i + (i + Ok)yh\n- 



We have 



£<!>' {nDbyh)rbDbV - — [ (j)' (V r,,yh)^ dx 

=eci>'{{i - oMt,^, + (1 + eMn){v',^+^ + 4, 



- —(I>{nyh\n+i) i nv' dx - —(p'{rbyh\n) [ nv'dx 
+ - ek)y'h\T,+, + (1 + ^fe)y;,|Tj - (1 - ^fe)0'(2y;,|T,+J - ^fe</''(2y;,|Tj]^^,+i k 



and 



e4>{rbDbyh) - — [ (piVr^yix)) dx 

=ls(^2m - 9Mn^, + ii + 9k)y'h\n) - {i + eM2y'^\n) - (i - 9k)cf>{2y'^\n^^)y (A.4) 

Subcase 3: If 6 = (4e, (4 + 2)£) , then rb = 2,bnTk= [4^, a;^] , 6 n Tk+i = [a;^, (4 + 2)e] , 
nDbV = v'^^^^ + v\^^^ and 
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We have 

ecf)' {nDhyh)nDbV - — [ (f)' {V r^yh)V r^v dx 
=£0'((2 - ek)y'f,\T,+, + eky'h\TM,+2 + v'e,+,) 
(t>'{nyh\Tk+i) I nv dx - —(t){nyh\Tk) ( nv'dx 

fb Jbf\Tk ^b Jbr\Tk+i 

=^1 [0'((2 - ek)y'n\n+, + ()ky'h\T,) - (i - eu)(t>' {2y'^\T,^^) - ek(t>' {2y'f,\TMu+i 
+ [<^'((2-^fe)y;,|T,+, +efeyklTj -</)'(2y;,|T,+J]^;^,+2}, 

and 

e(j){rbDbyh) - — I (j)iVrf,y{x)) dx 



n Jbnn 

1 



2 



s 2</.((2 - 9k)y'f,\T,+, + e^y'M - dkWh\n+,) - (2 - Okmy'hln) ■ (A.5) 



Case 2: b is across the left atomistic-continuum interface of an atomistic region. 
Subcase 1: If 6 = (4£,4+ie), then = 1, 6 n = (4^,3;^), b n = (x^, (4 + l)^), 
nDbV = v'i^_^_^ and 

nDbVh = (1 - Ok)yh\n+i + hy'hWk- 

We have 

ecj)' {rbDbyhjnDbV - |6 n fl3.\4>'irbDbnn^yh)Dbnn^v ~ ~ f 4>'{^nyh)^nV dx 



ct>'{{i - ek)y'n\n^, + Gky'h\n) - (i - Gk)4'' {y'h\n^,) - Okct>'{y'h\n 



vk+i, (A.6) 



and 

e4>{rbDbyh) -\br\ ^a.\rb(p{nDbnQ^y) - — [ <PC^rty{x)) dx 

^b Jbnnc 

=£(0((i - ek)y'h\n+, + hy'hlT,) - oMW,) - (i - 0kmyh\n+,)) (a.7) 

Subcase 2: If 6 = ((4 - 1)£, (4 -|- 1)£), then n = 2, b (1 Qc = ((4 - l)^,^;^), 6 n = 
(xj^, (4 + , rfeDfoV = v'^^^^ + and 

j-bAy/i = (1 - dk)yh\Tk+i + (1 + ^fe)y/ilTfc, nObna^yh = ^y^kk+i 

We have 

inDbyh)rbDbV -\bn iTbDbnQ.^yh)DbnQj^ ~ ~ I ^'{^nyh)^nV dx 

fb Jbncic 

=e| - ek)y',\n^, + (1 + 0fc)y;iTj - (1 - efc)0'(2y;jT,+J - 0k<P'{2y',\T,)]v'i^ 

+ [</.'((l - ek)y'n\n^, + (1 + ^fe)y;.|Tj - 0'(2y;,|rJ]^^;, - dk<l>' {y'hkM.+i, (A.8) 
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and 

£(f){rbDbyh) -\bn ^3\n(l){nDbnQ^y) - — [ (t){VrM^)) 

'''b Jbnnc 

=2e{m - ek)y'^\T,^, + (1 + 9M^,\n) - (1 - eM2y'^\T,^^) - (1 + 6M2y'^^\T,)). (A.9) 

Subcase 3: If 6 = (4e, (4 + 2)e), then rb = 2, bnnc = (4e,x^), 6 n f^a = (a^^, (4 + 2)e), 
VbDbV = v\^^^ + v\^^^, Dbm^v = 2^v'^^^2 + ^^^4+1, and 

nDbyh = y'h\Tk+2 + (i - ^k)y'h\Tk+i + ^ky'hWk nDbnn^yh = j:rf^y'h\Tk+2 + 2- el ^^l^^+i" 

We have 

e(p' {rbDbyh)rbDbV -\bn irbDbnn^yh)Dbnn^v - — [ (t>' {"^ nVhW nV dx 

+ {'^{^k\tu^. + (1 - ^fe)2/;.iT,+i + Qky'H\T,) - (1 - ^fc)0'(^y;.iT,+, + ^^^f^y;.iT,+j 

-0fe,^'(2y;,|Tj]^;,+i, (A.IO) 
Case 2: 6 is across the right atomistic-continuum interface of an atomistic region. 
Subcase 1: If 6 = (4£,4+i£^)) t^^en = 1, 6 fl ilc = (a^fe) (4 + l)^)) n fia = (4£)a;^), 
rfoDftV = v\^^^ and 

r-bAyfe = (1 - Ok)yh\n+i + hy'h\n- 

We have 



ecp' {rbDbyh)rbDbV - \b D Q^^lfp' {rbDbnn,yh)DbnQ^v ~ ~ [ (t>\^ nVh)^ r^v dx 



^'((1 - ek)y'H\n^, + ^kykln) - (i - ^fe)0'(y;.|T,+J - ekcl>' {ykln) 



4+1' 



(A.ll) 



and 

e(f){rbDbyh) - — I 4>C^ny{x)) dx 
fb Jbnnc 

=£0((i - ek)y'h\n+, + dkyhln) - [ Hv'hln) dx- [ </.(y;,|T,+i) dx 

Jbnn Jbnn+1 

=e{m - Ok)y'H\n^, + (^kykln) - 9Mh\n) - (i - Ok)Hy'h\n^,))- (A.12) 

Subcase 2: If 6 = (4e, (4 + 2)e), then rt = 2, 6 n = (aJ^, (4 + 2)£), 6 n f^a = (4£,a;^), 
rbL>fe?; = ?;;^+2 + ^4+i' nDbnn^v = v\^^-^ and 

''b-Dbyft = (2 - 6'fc)2/klTfe+i + 9ky'h\T^ and nDbna^yh = y'hln- (A.13) 

We have 

ecj)' {rbDbyh)rbDbV - |6 n ^3.\4>'irbDbna^yh)Dbna^v ~~i 4>'i^nyh)^rbV dx 

fb JbriCic 

=e!^ [</>'((2 - 0fc)2/;,|T,+, + Oky'^,\n) - Ok^P' {2y',^\T,) - (1 - 0k)<P' {2y'^\T,^^]v',^+, 

+ [0'((2 - eit)y;,|T,+i + eWhln) - ^'i'^y'h\n^M,+2, (a.m) 
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and 

e(l){nDbyh) - — I (/'(Vr.iy(x))dx 

'''b Jbnnc 

=e{m - 0k)y'h\n^, + Oky'hln) - oMhln) - (i - dkM^y.^J). (A.15) 

Subcase 3: If 6 = ((4 - l)e, (4 + l)e), then n = 2, bn^c = {x^, (4 + l)e), b n Q.^ = 
((4 - l)e,x^), nDbV = v'^^^^ + v'^^, nDbnn^v = j^v'^^^^ + j^v'^^ and 

nDbyh = (1 - &k)yh\Tk+i + hy'hln + y'hkk-i and nDbnn^yh = j^p^y/ilT^ + Y:^y'h\n-i- 

(A.16) 

We have 

Ecf)' {rbDbyh)rbDbV -\bf^ {nDbmJJh)Dbn^^v - — / (t>' nyh)^ nV dx 

JbnUc 

=£| [<^'((i - oMt,^, + e,y'^\T, + y'n\T,_,) - ^'{^/hln + i^/h\n.M, 
+ [0'((i - ek)y',\n^, + eky'„\T, + y'h\n_,) - Ok4>'{^^y'h\n + i^/h\n.,) 

-{l-ekWiy'^\T,^,)]v',^^„ (A.17) 

and 

e(f>{rbDbyh) - — [ H^rM^))^^ 
fb JbnQc 

=2£(0((i - 9k)y'h\n+, + ^ky'hln + y'hln.,) 

- (1 + (^'^^H^/hln + Y^/h\T,.,) - (1 - ^fc)0(2y;,|T,+J). (A.18) 

Appendix B. Approximation Properties 

In this section, we prove some approximation properties which we have used but are hardly 
found in standard text books. 

Lemma 11. Let v G C'^(M) n VF^'^(M) be a periodic function with [a,b] being one of its 
period. Let Vh be a Vi interpolation of v with respect to the nodes a < < xi < • • • < x„ < 
b < Xn+i = xq + (b — a) in [xo,Xn+i], subject to a constant, i.e., Vh{xk) = v{xk) + C, for 
A; = {1, 2, . . . , n + 1}, and is extended periodically with period b — a. Then the following estimate 
holds: 

II^ILf <ll^'llLf (B.l) 

[a,b\ [fiiO] 

where v' and v'^ denote the weak derivatives of v and respectively. 

Proof. First we note that, since v G C'^(M) and Vh is a Vi interpolation of v, the weak derivative 
of Vh on (xk,Xk+i) is defined by 

, _ v{xk+i) ~ v{xk) 
^h\^) — ■ 

Xk+l - Xk 

Since v G C''(M) is piecewise differentiable, we have 



v{xk+i) - v{xk) = / v'{t) dt, 
Jxk 
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where v' is the weak derivative of v. By the periodicity of and v' , and Cauchy-Schwarz 
Inequahty, we have 



"[a,6] 



n 



A;=0 ^'^"'"^ 
J a 



n 

1 

rb 

■dt 



'l|2 

a, 6] 



1^ \\L? 



Taking the square root on both sides gives the stated result. □ 



Lemma 12. Let v e C'^ {[a, b]) H W^'^{[a, b]) and IhV is the V\ function that interpolates v at 
the points a and b. We have the following inequality: 

\\v'-ihvy\\i. ,,<lb'lli^ • (B.2) 

(a,D) (a,o) 



Proof. Since v{a) = Ihv{a) and v{b) = Ihv{b), by the definition of I^v, we have 

j'b pb 

/ v' dx = / {Ihv)'dx, 

J a J a 

and equivalently, 

f\v' -{hvy)-idx = Q, 

J a 

where v' denotes the weak derivative of v on [a,h\. This shows that {Ihv)' is the best 
approximation of v' in the space of Vq functions as {Ihv)' is a constant. Therefore, by the 
property of best approximation, 

\\v'-{hvy\\l2 <\\v'-C\\l2 , (B.3) 

^ia,b) (a,b) 

for any constant C. In particular, if we choose C to be 0, the stated result holds. □ 
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Appendix C. Discrete Sobolev Inequalities on Non-uniform mesh 

In this section, we prove some discrete Soblev inequalities on non-uniform mesh that are used 
in the residual analysis for the external force. These results are extensions to the inequalities 
proved in [TJ Lemma A.l, Lemma A. 2, Theorem A. 4] on non uniform mesh. 

Lemma 13. Let g € M^, G and e^,ej > = 1, . . . , L, g' = (c/-)f=2 ^ M^-^ 

g'. := i = 2,...,L. If = 0, then 

L 



ffil < I^4l5fcl0i,fc, (C.l) 



h 

i=2 

where, h = Ya^^ 0i,fc = Y!1=i for k = 2, . . . ,i and (jji^k = YJt=k for k = i + 1, . . . , L. 
Proof. Let i G {1, . . . , L}, then 

L 

h\gi\ = \hgi - ^^'jOjl 

L L 

i=i i=i 

i-l L 
j=l j=i+l 



Since 



we have 



\9i-9j\ = I ^kd'kl, 
k=j+i 



i-l 



j=l k=j+l i=«+l k=i+l 

= E4i^^i(E4)+ E 4i.^i(E4) 

k=2 j=l k=i+l j=k 

L 

= Y^k\9'k\(pi,k- 

k=2 

Divide both sides by h, we obtain the stated result. □ 

Lemma 14. (Discrete Poincare's Inequality) Suppose that L > 1, s^,s^ £ M'^ with £^,ej > 0, 
\/i = 1,...,L. Let g e such that Yji=i^^9i = d' = i9'i)i=2 ^ ^^^'^ •'^c/i that 
g[ = Sj^-Sj—L_ Define Vq to be the set {1, . . . , L} and Vi to be the set {2, . . . , L}, then 

lL2max{maxi<j<Le°,max2<A:<Lel}^ii rr^ ^^ 

'"'"Vd.o) - 2 h '""'Vci^i)' ^ ^ 

for p G {1, oo}, where h = J2i=i ^i- 
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Proof. Using the result of Lemma 13, we have 











i=l k=2 



i=l 



i=l k=i+l 



L L 



k=2 i=l 



Since 



E < max e° ^ (l>i^k = ^max e° ^ 0i,fc + E 



i=l 



l<j<L 



fc-1 



1=1 



i=k 



and 



fc-i 



E '^^.'^ + E '^^.'^ < - 1) E + - - 1)) E 



1=1 



i=k 



l=k 



< [{k -l){L-{k- 1)) + {L-{k- l)){k - 1)] 



max e 

l<j<L 







< - max e°L^. 

2 l<j<L * 



Put these results together, we obtain the stated result for p = 1. For p = co, 

1 ^ 

\9i\ < v^4\9'k\(pi,k 



k=2 



< 



^4\9k\<Pi,k+ E ^k\9k\<Pi,k 



k=2 
L 



k=i+l 



k=2 



< 



1L2 



maxi<j<L e 







max eJgJ. 

2<k<L " " 



2 /i 

The stated result is obtained by taking the maximum of and over and 



□ 



Lemma 15. 

same as in Lemma 
fi = then 



(Discrete Friedrichs' Inequality) Suppose that L > 1, 
Let f £ such that f_ 



/l = 0, and f 



^ T>Q, T)2 are the 



P < - 

eOi-Do) 2 



(-^ ~ ^) o^™?r^ 1 max{e°,4}||/'||£!', 



(C.3) 



forp e {l,oo}. 
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Proof. For p = 1, 

L L-1 



For p = oo, 



and 



Thus we have 



i=l i=2 

L-1 



i=2 

= ^E^°OE(/.-/.-i)i + i 

i=2 j=2 j=i+l 

i=2 j=2 j=i+l 
i=2 j=l 



\fi\<^e]\f'j\ = {i-l) max^ e] max^ | /j | , 



J=2 
L 



\fi\ < > e]\f'j\ = (L — i) max ei[ max |/'|. 



max I f, | < mm(i — 1, L — i) max max I f'A 

id-Do 2<j<L ^ 2<i<L'^' 

< -(L — 1) max max If'l. 

- 2^ 2<j<L ^ 2<j<L' 



□ 



Remark 2. The bounds we have got here are not optimal as if e^'s and e^'s vary too much, 
taking the maximum of them in the inequahties could significantly reduce the sharpness of the 
estimate. However, for the analysis of this paper, such a bound is optimal enough to produce 
efficient error estimators and we leave the work of looking for optimal bounds to future work. □ 

Theorem 16. (hounds on the interpolation error) Let L > 1, e^, e^, G M^, with e^,e},£f > 
Vi = 1, . . . , L. Let f eR^ and F =G such that Fi = fi and 

i^. = /i + ^^^^(/L-/i) i = 2,...,L, (C.4) 

where h = Yli=2^i- define f = (/;)f^2 ^ M^-^ such that f[ = and f" = {f- )^~2 ^ 

M^~^ such that f" = — -, and F' and F" are defined in the same way. Let Dq, Di be the 
same sets defined in Lellmma 14 and V2 he the sei {2, . . . , L — 1}. Then, for p G {1, 00}, 
II, i-,|| lL^max2<j<L-ie-'max2<j<L-ie]max2<fc<L_ie| 

11/ - Fhyvo) < I 1 11/ hl,{V2y y^-^) 
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Proof. Let g = f — F, hy the definition of F, we have gi = gi = and 

L L L 



i=2 i=2 i=2 



By Lemma 15 



as = = 0, and by Lemma [T^ 



„ ,„ ^ 1 max{maxi<,<z. max2<fc<^-i ej}' „ 

Hi/ < 2 ii^' \m^2)^ 

as ^i=2^i9i ~ 0- Since -F" = 0, from which we know g" = f", the stated estimate holds. □ 
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